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Abstract 

We discuss how the main features of the recent LHC data on elastic scattering can be described by a 
QCD-inspired formahsm with a dynamical infrared mass scale. For this purpose new developments on a 
dynamical gluon mass approach are reported, with emphasis on a method to estimate uncertainty bounds in 
the predictions for the high-energy scattering observables. We investigate the effects due to the correlations 
among the fixed and free parameters involved and show that the band of predictions are consistent with the 
■ recent data from the TOTEM experiment, including the forward quantities and the differential cross section 
up to the dip position. 
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1. Introduction 

The TOTEM experiment has been designed to study elastic and diffractive scattering at the LHC, 
providing information on the pp total cross sections and the elastic differential cross sections at the energy 
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region 7-14 TeV and therefore expected to represent a crucial source of information for selecting models 
and pictures. First results at 7 TeV have already been published showing that no model predictions present 
full agreement with both total cross section and the differential cross section up to 2.5 GeV^ fllll- 

As a consequence reviewing processes are in order in the phenomenological context, and among the 
variety of approaches 0,13], the QCD-inspired models certainly play a fundamental role. However, although 
this class of models presents explicit connections with minijets and/or semihard concepts, where some kind of 
perturbative technique may be applied, the intrinsic soft character involved in the elastic scattering demands 
a quantitative connection with the non-perturbative QCD. 

With this aim (see also 0,0)13])' ^ novel interpretation for some fundamental parameters in eikonalized 
QCD-inspired models has been previously proposed, allowing a quantitative connection between elastic 
hadron scattering data and the nonperturbative regime [i, I3[l3]. This Dynamical Gluon Mass (DGM) 
approach is based on the possibility that the nonperturbative dynamics of QCD generates an effective gluon 
mass at the small momentum transfer region, which is strongly supported by recent QCD lattice simulations 



as well as by other phenomenological results |12l. Il3l. Il4|. In this context, the increase of the hadronic 



total cross sections is managed by the interplay between the gluon distribution function and a gluon mass 
scale and previous analysis on pp and pp elastic scattering has indicated good descriptions of the forward 
data up to = 1.8 TeV, including the differential cross section at small momentum transfer 

However, as commented above, the new results by the TOTEM collaboration demands further tests 
and developments in the phenomenological context and that is the point we are interested in here. As we 
shall show and discuss, even not using these new results as input information for quantitative fit procedures, 
general aspects of our formalism can be revisited and improved in the same energy region investigated before, 
namely between 10 GeV and 1.8 TeV (some preliminary ideas on these respect have already been discussed 



151]). Specifically, in this work we present new developments related to some formal and practical aspects 



of the approach with focus on: (1) quantitative account of the energy dependencies in the coupling constant 
and dynamical gluon mass; (2) fits based on improved data ensemble; (3) detailed analysis on the influence 
in the evaluated quantities associated with two fundamental parameters, the dynamical gluon mass scale 
and the soft Pomeron intercept; (4) estimation of bounds in all evaluated quantities associated with relevant 
physical intervals for the above parameters; (5) explicit reference to some calculational details. We shall 
show that these parameters play a crucial role in the energy region above that used in the fit procedures, 
which means they can not be fixed at ad hoc values without a fundamental physical justification. By 
considering relevant numerical intervals for the dynamical gluon mass scale and the soft Pomeron intercept, 
the question associated with the evaluation of the corresponding (statistical) uncertainty regions is addressed 
and a solution based on bounds for the high-energy predictions is proposed. Within these uncertainty bands 
we show that the recent data obtained by the TOTEM Collaboration at -^s = 7 TeV are reasonably well 
described, except the differential cross section beyond the dip position. Predictions at -y/s = 14 TeV and 
cosmic-ray energy region (Auger) are also presented. The advantages and some drawbacks associated with 
this formalism are discussed and further necessary developments outlined as well. 

The text is organized as follows. In Sect. 2 we treat the theoretical context, with explicit reference 
to novel formal account of the QCD scale and some calculational aspects. In Sect. 3 the analysis on the 
influence of the dynamical gluon mass and the soft Pomeron intercept are presented in detail, as well as the 
fit procedures, the proposed method for uncertainty inferences, the predictions at 7 TeV, 14 TeV (LHC), 
57 TeV (Auger) and a discussion on all the obtained results. The conclusions and some critical remarks are 
the contents of Sect. 4. 



2. Formalism 

2.1. Physical Quantities and Eikonal Representation 

Our analysis is based on three physical quantities expressed in terms of the elastic scattering amplitude 
yl(s,q), where q — \/—t, with s and t the Mandesltam variables: the differential cross section 



-^{s,q)=^\A{s,q)\'', 
2 



(1) 



the total cross section 



f7tot(s) =4rfmA(s,g = 0), (2) 

and the p parameter, 

Reyl(s,g = 0) 

From (1) and (2) the integrated elastic and inelastic cross sections can be obtained, 

da 



o-cl{s) ^ I ^{s,q)dq^, crin(s) = (Ttot(s) - crcl(s)- 
Jo dq^ 

In the eikonal representation (azimuthal symmetry assumed) the scattering amplitude is expressed by 

m 

A{s,q)=i bdbMqb){l-e'^^-''''^, (4) 
Jo 

where b is the impact parameter, Jq the zero-order Bessel function and x(s, b) the complex eikonal function 
in the impact parameter space. Eikonal models can be classified or distinguished according to different 
choices for the eikonal either in the impact parameter space or in the momentum transfer space, 

/"OO 

X{s,q)^ / bdbJo{qb)x{s,b), (5) 
Jq 

where it is connected with the specific dynamical assumptions that characterize each model. Under the 
physical condition Imx(s, b) > this representation constitutes an automatic unitarized framework suitable 
for theoretical developments. In special, in the case of QCD based or inspired models the eikonal in the 
momentum transfer space is expressed in terms of hadronic form factors and elementary cross sections and 
the main point is how to connect these quantities with parton-parton scattering processes. Since in these 
models the rise of the total cross section at high energies is driven by gluon-gluon semihard scattering, this 
is the main process to be investigated. 

2.2. Dynamical Gluon Mass and Gluon-Gluon Cross Section 



In addition to lattice QCD the Schwinger- Dyson equations constitutes a fundamental nonperturba- 
tive formalism in the investigation of hadron physics, specially in what concerns confinement and dynamical 
chiral symmetry breaking. Infrared finite solutions for these equations have been obtained by Cornwall by 
means of the pinch technique leading to a gluon propagator which carries a dynamical mass 
review on this dynamical gluon mass generation is presented in [23|, discussions and more details can be 
found in Q and references therein. Here we display the main formulas connected with the nonperturbative 
gluon-gluon cross section, stressing a formal improvement regarding the previous work [9]. 

In the DGM approach the nonperturbative integrated elastic cross section for the scattering gg — >■ gg, 
expressed in terms of the subprocess energy s, is given by Q 



g g g ^ Bin ' ^ 



4M2s[s-M2]2 I M2 



(6) 



where = as{s) and Mg = Mg{s) are the running coupling constant and the dynamical gluon mass, 
respectively. From the Cornwall solution for the gluon propagator in the case of pure gauge QCD, these 



quantities are expressed by [17|, |18|, |1 



47r 

a (s) = — (7) 

' /3oln[(s + 4M2(s))/A2]' 

3 
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-12/11 



(8) 



where /3o = 11 — |n/ (n/ is the number of flavors), A = Aqcd and rrig = Mg{0) is the gluon mass scale. A 
fundamental consequence of this result is its extension to perturbative QCD. In fact, in the limit s ^ A^ 
we have Mg{s) — >■ and matches the one-loop perturbative QCD one 

In this context, by considering nf — A and A = 284 MeV, the only unknown quantity in the evaluation 
of o-gg(s) is the gluon mass scale rUg, the IR value for Mg js). Phenomenological analyses suggest typical 



values of TTin in the range from 300 MeV up to 700 MeV 0, H El 0, 0, bj, UJ, Il4f and in |9|| it has been 



inferred rUg = 400 MeV. To investigate the relevant numerical interval for uig and the effects in the analysis 
of elastic hadron scattering is one of the main points we are interested in here. 

It should be noted that in previous works d, 0, [13, 13 1 the dependence of Mg on the energy of the 
subprocess gg — >■ gg has been neglected as a first approximation, that is, in Eq. ^ the IR value has 
been used in place of Mg{s). Therefore, the integrated cross section as given by Eq. (6) generalizes the one 
considered in 0, [l^l , enlarging the physical meaning of the dynamical gluon mass and the efficiency of 
the formalism in the data analyses, as shown in Sect. 3. Another parametrization for Mg{s) is discussed in 
Appendix A. 



2.3. Eikonalized Approach 

At high energies the soft and the semihard components of the scattering amplitude are closely related 
23I 24 1, and it becomes important to distinguish between semihard gluons, which participate in hard parton- 
parton scattering, and soft gluons, emitted in any given parton-parton QCD radiation process. Hence a 
formalism based on QCD has to incorporate soft and semihard processes in the treatment of high-energy 
hadronic interactions and, more importantly, has to bring up information about the infrared properties of 
QCD. 

In the phenomenological context, the developments of the mini-jet or semi-hard QCD picture Q con- 
stitute an interesting step in the search for connections between hadron scattering and QCD subprocesses, 
associated with gluon-gluon [gg), qua rk-g luon {qg) and quark-quark {qq) interactions. That is the case in 
the formalisms by Margolis et al. [25|, |26|, [27 1 and the more recent version sometimes referred to as Aspen 
model, by Block et al. }^^^. However, despite its efficiency in the description of the forward quantities in 
the elastic channel, as well as differential cross sections at small values of the momentum transfer a serious 
drawback with the Aspen model concerns the absence of the expected connections with nonperturbative 
QCD. Moreover, and as a consequence, two fundamental quantities, the infrared mass scale, denoted mo, 
and the coupling constant a^, are unknown parameters fixed to ad hoc values 600 MeV and 0.5 respectively, 
in order to obtain best fits in data analyses 2^ 2^. 

In this respect, the main ingredient in the DGM approach [^, Q has been the novel physically moti- 
vated interpretation for both quantities through the gluon mass scale as a natural regulator of the infrared 
divergences associated with the semihard gluon-gluon subprocess cross sections, leading therefore to direct 
connections with nonperturbative QCD. The general formulation is largely based on the eikonal structure 



of the Aspen model [28|, |29| , but with the above well founded connection, as well as differences regarding 



free/fixed parameters, constraints and specific parametrization. The point is to consider Eqs. (6 - 8) for 
the gg subprocess and the IR values as{s = 0) and M^(s = 0) = to^ in substitution to as and toq in aU 
elementary subprocesses. 

Specifically, for pp and pp elastic scattering, the complex eikonal in the impact parameter space is 
expressed in terms of even (e) and odd (o) functions. 



Xpp{s, b) = Xe{s, b) + Xo{s, b), Xpp{s, b) = Xe{s, b) ~ Xo{s, b). 



(9) 



The odd part is expected to account for the difference between both reactions at low energies and to 
vanish in the high-energy limit. It is therefore parametrized as a Reggeon contribution and in a factorized 
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form, with the inclusion of the dynamical gluon mass ingredients and the odd prescription, 



Xois,b)^kCo^e'''^^W{b;fio), (10) 

where 

(11) 

W{b; fj,o) is the Fourier transform of the dipole form factor 

Wib;fio)^^M'K3M (12) 
9D7r 

(K^ is the modified Bessel function of second kind), and fio and Co are free fit parameters. 

The even part is supposed to be expressed in terms of contributions from elementary interactions, namely 
gluon-gluon, quark-giuon and quark-quark, in a factorized form at the constituent level, 

Xe{b,s) = Xqq{b,s)+Xqg{b,s)+Xgg(b,s) 

= i[aqq{s)W{b;fiqq) + aqg{s)W{b;tiqg) + agg{s)W{b;figg)], (13) 

where, for ij = qq,qg,gg, W{b,fiij) are overlap functions with the same structure of Eq. (12), /ijj are free 
fit parameters and aij are the corresponding elementary cross sections. 

For the sub-processes involving quarks, the parametrizations for the cross sections follows from approxi- 
mate forms of distribution functions involving quarks and gluons at small x region, with the even prescription 
in the corresponding eikonal: 

^,,(.s) = fcC,,^e-/4, (14) 



aqg{s) ^ k I Cgg + C' 



qg 



In 



(15) 



where Cgg, Cqg, Cgg, are free fit parameters. The convenience of the general structure adopted in the 
parametrizations (fTU|) . ([T^. and ((15]) can be traced back to the early work of Margolis et al. 25, 2^ 2^ 
and references therein. 

The full gg cross section is expressed by 



99 



4»Tl2/s 



dTFggir) aggis) 



where r — s/s, Fgg{T) is the convoluted gluon distribution function, 

1 dx 



f dx / T 

Fggir) ^ [g ® g]{T) = — g{x) g [- 



(16) 



(17) 



and agg{s) represents the subprocess gg — )■ gg cross section, given here by the nonperturbative result (6). 
For the gluon distribution function, since in the DGM approach the small- a; semihard gluons play a central 
role, the phcnomcnological parametrization used in [28] is also considered here, 



{l~xf 



(18) 



which simulates the effect of scaling violations in the small x. Here J can be associated with the soft Pomcron 
intercept ( J = 1 + e). Notice that the ansatz (|18p satisfies the normalization condition dxxg{x) — A, 



where A is a constant [3J|. If A = 1/2, Ng = (6 — e)(5 — e)...(l — e)/240. The power 5 appearing in the term 
(1 — a;) is suggested by dimensional counting rules |34l. l35l. l36l|. 

The integration of Eq. (16), taking account of the analyticity properties, demands some comments. The 
odd contribution to the eikonal, Eq. (10) and those from qg and qq subprocesses to the even eikonal, Eqs. 
(14) and (15), are closed analytical forms, so that the corresponding prescriptions are straightforward to be 
applied, leading to the generation of the real and imaginary parts. On the other hand, the gg contribution, 
Eq. (16), cannot be analytically evaluated and any numerical approach, embodying analyticity, demands 
some assumptions. We have used two independent methods. In the first one generated numerical points for 
the real primitive of the integrand in (16) have been parametrized by polynomials in Ins and then the even 
prescription has been applied. In the second one we have used first order derivative dispersion relation for 



even functions [37[ applied as 



Recrgg(s) w -^-j^^^^'^ggi^)^ (19) 



and numerical derivation of the generated points. In (19), Imcrgg(s) corresponds to the generated numerical 
points for the real primitive of (16). Both methods lead to the same results, except for negligible differences 
at low energies. 

With all the above ingredients the formalism is complete. Specifically, it has two fundamental parameters, 
the dynamical gluon mass scale nig and the soft Pomeron intercept e. Once mg and e are chosen, the model 
has 9 free parameters to be determined by fits to the experimental data: 5 normalization constants associated 
with the odd contribution and elementary cross section (Co, Cgg, Cqg, C'^g, Cgg) and 4 coming from the 
corresponding form factors (/io, /igg, /igg, /^gg). That should be contrasted with 6 fixed and 6 fit parameters 
in the Aspen model [Ulli]. 

In what concerns the momentum transfer space, the essential ingredient in the structures of the eikonals 
(9), (10) and (13) is the combination of four dipole form factors, of type (12) in the impact parameter space 
(in fact, a linear combination at each fixed energy). Therefore, the dependence of the physical quantities on 
the momentum transfer is a direct consequence of these choices and in this respect the following comments 
are in order. Historically, the dipole ansatz may be traced back to the conjecture by Wu and Yang that the 
distribution of hadronic matter inside a hadron is similar in electromagnetic and strong interactions (soj . 
Specifically they conjectured that the elastic pp differential cross section might be proportional to the fourth 
power of the proton charge form factor, that is the form factor measured in electron-proton scattering. In 
the context of the Chou-Yang model [3l|, with an input for the form factor, the differential cross section 
can be evaluated. The use of a parametrization as a sum of exponential in q^ by Chou and Yang (3^ and 
the dipole electric form factor by Durand and Lipes [33] . 

+ M^ = 0.71GeV^ 

led to the prediction of the diffractive pattern in pp differential cross section and, more importantly, the 
correct position of the dip, as experimentally observed latter. However, the dipole parametrization predicts 
multiple dips and peaks that have never been observed in the experimental data. In our case (as well as 
in the Aspen model), the oscillatory behaviour, characteristic of the dipole structure can be attenuated by 
the combination of four different dipoles, at least up to the intermediate region of the momentum transfer 
(dip-bump structure). However, beyond this region, changes of slopes are expected, as indicated by the 
results presented in the following section. This effect seems not to be present in the published plots by the 
TOTEM Collaboration [l[, suggesting that other form factors should be regarded in order to describe the 
high-t data. We shall return to this point in the conclusions and final remarks. 

3. Testing parameters: Fits and Results 

Despite the suitable formal connection with nonperturbative QCD and the efficiency in the description of 
the experimental data, one important aspect to be investigated in the DGM approach concerns the influence 
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of the fundamental parameters in the evahiated quantities (fits and predictions): nig and e. They could 
certainly be fixed, as done for example with the corresponding parameters (and others) in the Aspen model 
[28| . However, as we shall show, relevant numerical intervals for rUg and e affect the model predictions at 
the highest energies and that puts serious limits on the reliability of models with ad hoc values for the fixed 
parameters involved. Our goal here is to point out the relevance of this aspect and the fact that it demands 
some solution or a clear explanation, prior to consider these models as physically well founded. To this 
end we first discuss some aspects of the data ensemble used here (which differs from the one considered 
in [1, 0) and outline the fit procedure (Subsection 3.1). Then we treat in some detail the selection of the 
relevant intervals for the parameters and the proposed methodology for uncertainty inferences (Subsection 
3.2), followed by a critical discussion on the obtained results (Subsection 3.3). 



3.1. Data Ensemble and Fit Procedure 

In the previous work Q the data ensemble used in the fits includes only forward quantities, utot^ P 
parameter and the slope parameter B up to 1.8 TeV. Here we shall consider a different ensemble for the 
reasons that follows. 

Even in the very forward direction the evaluation of the slope parameter depends somehow on the 
momentum transfer interval considered (at least two points), which renders difficult comparison among 
results from different experiments with different intervals. That is clearly illustrated in the comparison of 
the recent result by the TOTEM Collaboration with the existing data [38| that includes the pp2pp result 
at 200 GeV [1^. Moreover, in a model context the slope is evaluated at one point, = 0, and therefore 
some bias may always be present in the calculation. For that reason, we shall not consider here the forward 
slope as input quantity in the data reductions. 

On the other hand, as constructed, the formalism is expected to be applicable in the soft region, which 
in terms of the differential cross section means intervals in momentum transfer roughly up to 1 GeV^, or 
even beyond. The presence of four dipole form factors in the eikonal (Subsection 2.3) corroborates this 
expectation. Hence we include in the ensemble to be fitted the differential cross section data available at the 
highest energies, namely pp scattering at 546 GeV and 1.8 TeV, covering the region up to ~ 1.5 GeV^ 
and ~ 0.6 GeV^, respectively. 

Since the numerical tables on the differential cross section at 7 TeV are not yet available and a second 
measurement of the total cross section (through a luminosity independent method) is expected, we shall not 
consider here the recent TOTEM results as input in the fits. However, the numerical information available 
(and published) will be referred to as a guidance in the discussion of our predictions. 

Our ensemble therefore includes the pp and pp data on atot and p above 10 GeV up to 1.8 TeV, together 
with t he p p differential cross section data at 546 GeV and 1.8 TeV. The data have been collected from the 
PDG [40| and Durham 41 1 sites and in all cases the statistic and systematic errors have been added in 
quadrature. The data reductions consist in assuming some input values for m^ and e and through Eqs. (1) 
to (3) to fit the above data ensemble. For the fit program we have used the class TMinuit of the CERN 
ROOT Framework [i^l and the MIGRAD minimizer, setting the confidence level at 90%. 



3.2. Testing Parameters 

Preliminary results with the above ensemble and fit procedure have already been presented in our previous 
analysis [15| . The main conclusion has been that by considering narrow intervals for rUg and e all the fitted 



data are quite well described in all cases. However, at higher energies the predictions are sensitive to each 
numerical interval considered. Here, we address the question on how to infer quantitative uncertainties in 
the model predictions taking into account relevant intervals for rrig and e. It is important to stress that our 
point does not concern standard propagation from the errors in the fit parameters, as given by the error 
matrix (variances and covariances). Our focus is only on the effects of physically motivated intervals for rUg 
and e. 

However that is not an easy task, due to the numerical integrations involved, the absence in this case 
of a statistical method to extract uncertainties and the strong correlations among m^, e and all the free fit 
parameters, as we shall shown. What we propose here concerns a method to infer bounds in the evaluated 
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quantities in accordance with relevant intervals for mg and e, in a similar way as has been recently done by 
Achilli et al. in the context of minijet models [i^ . 

Due to the correlations among e, uig and the fit parameters, a relevant interval for one of the fundamental 
parameters depends on the relevant interval for the other one. For that reason we first treat e and then irig. 
For the soft Pomeron intercept, as in [l^, we consider the results from the detailed analyses on the extrema 



and constrained bounds for the intercept, since they are based on both scattering data [4j| and scattering 
together with spectroscopy data [i^, respectively. With these bounds, namely eiower = 0.080 and eupper = 



0.090 [4J,|45|, we consider as relevant interval all real values in this range, including the bounds. Thus we 



use the following notation for intervals with that meaning (equally likely values): 

e : [0.080, 0.090]. 

For the dynamical gluon mass scale, as commented in Subsection 2.2, phenomenological analyses suggest 
Trig values in the range from 300 MeV up to 700 MeV. Taking into account the interval for e, as we shall 
show, the following interval for m^, larger than that considered in [l5| . and with the above notation and 
meaning, is adequate to our purposes: 

mg(MeV) : [300, 600]. 

Once selected the relevant intervals, the next step is to investigate their effects in the fit results and 
evaluated quantities, looking for the corresponding uncertainties in these quantities. At this point it is 
important to stress that, in principle, any numerical pair (e, rrig) in the selected intervals is equally likely 
and therefore we can not attribute central values with uncertainty, as for example 0.085 ± 0.005 for e and 
450 ± 150 for mg. That is, there is no meaning in assuming that a most probable value is any closer to the 
middle of this range that it is to either endpoint. Moreover, it is also not possible to test the effects of all 
real pairs (e, rUg). For those reasons we have considered a method consisting of the following steps: 

1. From different tests with real pairs (e, rUg) we have selected relevant values for each parameters in 
each interval, with the result: 

e : 0.080, 0.085, 0.090 
m„ : 300, 400, 500, 600 MeV 



2. For each of the above pairs (e, rUg), therefore 12 variants, we have performed the data reductions 
(Subsection 3.1), determining in each case the per degrees of freedom (DOF) as a test for goodness 
of fit. 

3. Based on the above statistical information and in the global description of the fitted data, we have 
determined an optimum value for one of the parameters together with the relevant interval for the 
other one. 

4. By fixing one parameter in its optimum value the extrema values of the interval corresponding to the 
other parameter have been used as upper and lower bounds in the evaluated quantities. 

5. The band limited by the above upper and lower bounds has been considered as an estimation of the 
uncertainty region in the evaluated quantities. 

With this strategy we do not take account of the simultaneous influences of both intervals, but only the 
effect of the interval in one parameter for a fixed value of the other parameter. We shall return to this point 
in subsection C. 

The values of the x^/DOF for each one of the 12 variants are displayed in Table I and plotted in Fig. 
1, in terms of the dynamical gluon mass scale and for each value of the soft Pomeron intercept. Since we 
have 318 DOF the best results demand x^/DOF closest and/or below 1. Figure 1 plays a central role in our 
analysis since we can extract the following main conclusions: 
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Table 1: Values of the x^/DOF for 318 DOF obtained in all fits. 







rUg 


MeV 




e 


300 


400 


500 


600 


0.080 


1.160 


0.9487 


0.9690 


0.9984 


0.085 


0.9899 


0.9536 


0.9536 


1.432 


0.090 


1.001 


0.9597 


1.070 


4.339 



1. For nig = 600 MeV the cases e = 0.085 and 0.090 are excluded on statistical grounds. 

2. For e = 0.080, despite the reasonable result at nig = 300 MeV, the whole interval in nig corresponds 
to quite good statistical results. 

3. For nig = 400 MeV the whole interval for e corresponds to quite good statistical results. 

4. For both parameters considered simultaneously the intervals e: [0.080, 0.090] and nig (MeV): [400, 

500] correspond to acceptable statistical results. 

Despite the last conclusion, in order to develop here a feasible and economic method for uncertainty 
inferences, we do not consider both intervals simultaneously. For individual intervals, based on the first 
three statistical conclusions and after checking the descriptions of the fitted data in all cases analyzed, we 
propose the following optimum solutions, with the interval notation and meaning defined before: 

- For fixed e = 0.080 relevant interval nig (MeV): [300, 600]; 

- For fixed nig = 400 MeV relevant interval e: [0.080, 0.090]. 

We note that the asymmetrical character of this solution is a consequence of the statistical analysis we 
have considered. Moreover, as we show in what follows, this solution allows us to estimate the uncertainties 
in the evaluated quantities in a consistent way and that is the point we are interested in here. 

Table 2: Values of the fit parameters for e = 0.080. Co, Cqq, Cqg, C'^g and Cgg are dimensionless and 
GeV. 



mg (MeV) : 


300 


400 


500 


600 


Co 


1.170±0.026 


3.03±0.40 


5.61±0.24 


10.12±0.41 


Cqq 


3.8498±0.0036 


10.7±1.4 


30.65±0.72 


21.2±1.9 


CqgixlO-^) 


1.8993±0.0082 


8.74±0.59 


13.414 ±0.034 


42.70±0.62 


C;,(xl0-2) 
C,,(xl0-3) 


1.8174±0.0068 


4.51±0.62 


18.842±0.060 


18.55±0.74 


0.979±0.019 


3.79±0.17 


6.832±0.057 


18.20±0.26 


Mo 


0.1862±0.0041 


0.41±0.17 


0.196±0.021 


0.622±0.053 


l^gg 


0.6463 ±0.0055 


0.651±0.066 


0.6314±0.0031 


0.6427±0.0020 


0.156218±0.0023 


1.32±0.16 


1.574±0.047 


1.03± 0.15 


^Iqg 


0.8368±0.0027 


0.838±0.044 


0.8214±0.0020 


0.8377±0.0018 



Let us first consider the case e = 0.080 and nig (MeV): [300, 600]. The values of the fit parameters for 
nig = 300, 400, 500 and 600 MeV are shown in Table II. With the extrema values nig = 300 MeV and nig 
= 600 MeV for each evaluated quantity we plot the upper and the lower curves delimiting the region of 
uncertainty represented by a band, in the figures that follows. The results for atot{s), p{s) are shown in Fig. 
2 and those for the pp differential cross sections at 546 GeV and 1.8 TeV in Fig. 3. 
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Figure 1: Reduced chi-square in terms of the dynamical gluon mass scale for different values of the soft Pomeron intercept. 
The straight line segments are drawn only to guide the eyes. 



For further discussion, we have included in these figures some data that did not take part in the data 
reductions. In the plot of atot (Fig- 2. a) we have included the recent result by the TOTEM Collaboration 
at 7 TeV and in the differential cross section at 1.80 TeV (Fig. 3.b) the data at 1.96 TeV 

The predictions for the pp differential cross section data at 7 TeV and 14 TeV are displayed in Fig. 4. In 
the former case (7 TeV) the straight line segments and the point indicate recent experimental results obtained 
by the TOTEM Collaboration [l|, Q (see Table V) : the optical point with the slope at the diffraction peak 
(dotted segment), the dip position (dashed vertical segment) and the value of the differential cross section 
at = 0.7 GeV^ (black square). 

For the case rUg = 400 MeV and e: [0.080, 0.090], the values of the fit parameters for e= 0.080, 0.085 
and 0.090 are displayed in Table III. Here, with the extrema values e — 0.080 and e — 0.090 we plot the 
upper and the lower curves delimiting the region of uncertainty represented by a band in the figures that 
follows. The corresponding results for crtot(s), p{s) and the differential cross sections are shown in Figs. 5, 
6 and 7. 

3.3. Discussion 

The main point in this section has been to propose estimations of the uncertainties in the evaluated 
quantities associated with relevant intervals for the fundamental parameters e and nig. As commented before 
we did not address the uncertainties associated with both intervals simultaneously but only an individual 
interval for a fixed parameter. Although that may represent an underestimation of the effective uncertainty 
it certainly gives more information on the reliability of the results than what could be obtained by fixing 
both parameters (or 6 as in (28j'l. In addition, in practice one band is practically contained in the other, so 
that the larger bands represent a reasonable result for our purposes. 

The results presented in Tables I - III, figures 2, 3, 4 and 5, 6, 7 lead to some qualitative conclusions: 
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Figure 2: Total cross sections and p parameter for fixed e = 0.080, upper bound for mg = 600 MeV and lower bound for mg 
= 300 MeV. 




0.5 1 1.5 2 ^"'o 0.5 1 1,5 2 

Itl(GeV') Itl(GeV') 



Figure 3: Differential cross sections at (a) 546 GeV and (b) 1.8 TeV (togctfier with the 1.96 TeV data) for fixed e = 0.080, 
upper bound for mg = 600 MeV and lower bound for rrig = 300 MeV. 
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Figure 4: Predictions for the differential cross section at 7 TeV (a) and 14 TeV (b) for fixed e = 0.080, upper bound for rrig 
= 600 MeV and lower bound for rUg = 300 MeV. At 7 TeV the dotted segment indicate the optical point and slope at the 
diffraction peak, the dashed vertical segment the dip position and the black square the differential cross section at = 0.7 
GeV^, as measured by TOTEM Collaboration (see Table V). 




Figure 5: Total cross section and p parameter for fixed rrig = 400 MeV, upper bound for e = 0.090 and lower bound for e = 
0.080. 
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\\\m 



Figure 6: Differential cross sections at (a) 546 GeV and (b) 1.8 TeV (together with the data at 1.96 TeV) for fixed nig = 400 
MeV, upper bound for e = 0.090 and lower bound for e = 0.080. 
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Table 3: Values of the fit parameters for rrig = 400 MeV. Co, Cqq, C'qg, C'gg and Cgg arc dimensionless and /lo, /J-gg, fJ-gg, fJ-gg 
are in GeV. 



0.080 0.085 0.090 



Co 


3.03±0.40 


3.10±0.48 


3.11±0.42 


Cqq 


10.7±1.4 


10.5±1.2 


10.2±1.1 


CqgixlO-^) 


8.74±0.59 


8.63±0.47 


8.66±0.41 


C;,(xlO-2) 
(X 10-3) 


4.51±0.62 


4.68±0.50 


4.69±0.45 


3.79±0.17 


3.62±0.14 


3.49±0.12 




0.41±0.17 


0.44±0.18 


4.45±0.16 




0.651±0.066 


0.6500±0.0065 


0.6496±0.0063 


Hqq 


1.32±0.16 


1.30±0.14 


1.27±0.16 




0.838±0.044 


0.8367±0.0040 


0.8367±0.0037 



i. In addition to the good quality of the statistical results in terms of x^/DOF (Table I), Figs. 2, 3 and 
4 show that the visual description of all reduced data is quite good. In this energy region (up to 1.8 
TeV) the inferred error bands are not significant, which means that for any fixed values of e and nig 
considered in the analysis, all data are quite well described. The correct reproduction of the high- 
energy differential cross section data up to ^ 1.5 GeV^ (Figs. 3 and 6) seems to be a remarkable 
result in the context of QCD-inspired models. 

ii. At higher energies (above those considered in the data reductions) the effects of the relevant intervals 
are significant since the general trend is an increase in the estimated error band above 1.8 TeV (Figs. 
2 and 5). 

iii. The error bands for fixed rrig (400 MeV) are narrower than those for fixed e (0.080) and the former 
are, in general, contained in the latter. That is a consequence of the difference in the relative intervals 
involved, ^ 6 % for e and ^ 30 % for nig, as well as the role played by each parameter in the formalism 
and data reductions. 

iv. Table II and III illustrate the correlations among e, nig and the fit parameters. We see that this 
correlation is stronger in the case of e fixed than in the case of nig fixed, an effect also associated 
with the difference in the relative intervals. Moreover, note that in both cases the correlations are 
stronger among the parameters associated with the normalization constants (the C's) than among 
those coming from the form factors (the /i's). 

V. Taking into account the error bands, the total cross section result by the TOTEM Collaboration at 7 
TeV [4] is quite well described in both cases (Fig. 2. a and Fig. 5. a) and a reasonable consistency can 
be observed in the differential cross section up to the dip position (Fig. 4. a and Fig. 7. a). 

vi. For the differential cross section at 14 TeV a change of curvature is predicted around q^ ~ 2 GeV^ in 
both cases (Fig. 4.b and Fig. 7.b). 

Let us now present and discuss some quantitative results and predictions with the 12 variants considered, 
as well as the inferred upper and lower bounds for the evaluated quantities. For each one of the 12 variants 
the results for atot, P, (^in and a ei jo tot at 7, 14 and 57 TeV are displayed in Table IV. As commented before, 
with the extrema values of e or nig in each relevant interval, we evaluate the upper and lower bounds for all 
evaluated quantities. Table IV shows that in fact, for the values of e and nig inside the relevant intervals 
the numerical values of the evaluated quantities lies inside the upper and lower bounds. Therefore, in what 
follows we shall give our predictions in terms of the band corresponding to these upper and lower bounds 
in each case (e or nig fixed). Since these bands have the same meaning as the relevant intervals (equally 
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Table 4: Predictions for pp integrated cross sections and p parameter at the LHC and Auger energies for each tested value of 
rrig and e. For each physical quantity the extrema values define the uncertainty band. The cross sections are in mb. 



physical 


e = 0.080 










TUg = 400 MeV 






quantity 


mg (MeV) : 


300 


400 


500 


600 


e: 0.080 


0.085 


0.090 


(Tin (7 TeV) 




70.0 


72.3 


70.6 


74.6 


72.3 


72.6 


73.0 


a„ (14 TcV) 




76.7 


79.8 


77.5 


82.9 


79.8 


80.3 


81.0 


a™ (57 TeV) 




91.4 


96.1 


92.8 


100.9 


96.1 


97.2 


98.4 


<Jtot (7 TeV) 




93.2 


96.9 


94.0 


100.2 


96.9 


97.4 


98.0 


atot (14 TeV) 




103.9 


108.8 


104.9 


113.4 


108.8 


109.6 


110.6 


atot (57 TeV) 




127.8 


135.6 


129.5 


143.1 


135.6 


137.3 


139.4 


fTeiMot(7TeV) 




0.2489 


0.2539 


0.2489 


0.2554 


0.2539 


0.2546 


0.2551 


a,i/atot{U TeV) 




0.2618 


0.2665 


0.2612 


0.2690 


0.2665 


0.2673 


0.2676 


a'ei/o'tot(57 TeV 




0.2848 


0.2834 


0.2834 


0.2949 


0.2834 


0.2921 


0.2941 


p{7 TeV) 




0.1225 


0.1321 


0.1245 


0.1409 


0.1321 


0.1346 


0.1376 


p{U TeV) 




0.1189 


0.1272 


0.1208 


0.1349 


0.1272 


0.1299 


0.1330 



likely values), we shall use the same notation, namely the bounds inside brackets and separated by a comma: 
[lower-bound, upper-bound] . 

Let us consider first the pp elastic scattering at 7 TeV with focus on the published results by the TOTEM 
Collaboration [ll,0|- Our predictions (bands) for several physical quantities are displayed in Table V, together 
with the corresponding results by the TOTEM Collaboration. From that Table and within the error bands 
we are led to the conclusions that follows. 

• For e — 0.080 the predictions are: (a) consistent with all cross sections, optical point and dip position; 
(b) barely consistent with the forward slope; (c) not consistent with the slope near the dip position, the 
differential cross section at — 0.7 GeV^ and the exponent n in the power law at large momentum transfer. 

• For TTig = 400 MeV the predictions are: (a) consistent with all the cross sections and the optical point; 
(b) barely consistent with the dip position; (c) not consistent with both slopes, the differential cross section 
at = 0.7 GeV^ and the exponent n. 

How good or bad are these results? In [ij the TOTEM Collaboration displays in Table 4 the predictions 
from some representative phenomenological models compiled in [47] . Comparison of our results with these 
predictions shows that despite all the above inconsistencies, our efficiency in the description of the experi- 
mental data is compatible with that presented by all quoted models. It is also interesting to note in that 
Table that our wrong prediction for the exponent in the power law at large q^, namely n ~ 10.5, is in plenty 
agreement with the prediction by Block et al.: n — 10.4. This agreement may be related with the choice 
of four dipolc form factors, present in both models. Moreover, the first TOTEM data on differential elastic 
cross section at 7 TeV has been measured in specialized runs of LHC and in the momentum transfer range 
of 0.36 < —t < 2.5 GeV^. Notice that none of the representative models for elastic scattering can fully 
reproduce its overall features; in special, the dip position, dip depth and large-t behaviour [IJ. Particularly 
in the range 1.5 < —t < 2.0 GeV^ the data presented a power law pattern proportional to ~ which is 

in poorly agreement with general model predictions. Interestingly enough, this specific power law behaviour 
was first envisaged by Donnachie and Landshoff 48] in a model which accounts for the contribution of three 
gluon exchange in elastic scatte ring at high energies. Yet, in the QCD context, the constituent interchange 
model by Lepage and Brodsky 49] predicts a different power law behaviour, namely ~ The DGM 

model results at the large-i region are consistent with the latter. We shall return to this point in our 
conclusions. 

At other energies (14 and 57 TeV) the bands in our predictions also correspond to the extrema values 
for each physical quantity, as given in Table IV, for the cases e = 0.080 and rUg = 400 MeV. For example. 
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at 57 TeV we predict atot (mb): [127.8, 143.1] for e = 0.080 and atot (mb): [127.8, 139.4] for mg = 400 
MeV. We note that these results are consistent with recent estimation of this quantity through an analytical 
parametrization, atot = 133.4 ±1.6 mb [Hol ]. 



Table 5: TOTEM results at 7 TeV with systematic and statistical error added in quadrature and the uncertainty band 
predictions. 



Physical 


TOTEM 


e = 0.080 


rUg = 400 MeV 


quantity 


results 


rUg (MeV): [300, 600] 


e: [0.080, 0.090] 


B(0.36 < t <0.47 GeV^) (GeV"^) [1] 


23.60 ± 0.64 


[19.8, 22.4] 


[20.8, 21.6] 


\tMp\ (GcV) \l] 


0.53 ± 0.01 


[0.51, 0.54] 


[0.51, 0.52] 


n in \t\~" (1.5 < t <2.0 GcV^) [1] 


7.80 ± 0.32 


[10.1, 10.6] 


[10.5, 10.5] 


f^(|t| = 0.7) (mbGeV-^) [1] 


2.70^[;i^xl0-^ 


[3.0, 4.3] X 10-2 


[3.8, 4.1] X 10-2 


del (mb) \2] 


24.8 ± 1.2 


[23.2, 25.6] 


[24.6, 25.0] 


(Tin (mb) [2] 




[70.0, 74.6] 


[72.3, 73.0] 


atot (mb) [2] 


98.3 ± 2.8 


[93.2, 100.2] 


[96.9, 98.0] 


B(0.02 < t <0.33 GeV^) (GeV^^) [2] 


20.10 ± 0.36 


[18.6, 19.8] 


[19.3, 19.4] 


^^\t=o (mbGcV-2) \2] 


504 ± 27 


[451, 523] 


[488, 500] 



At last, let us discuss our two proposed independent solutions for the inference of uncertainties. Guided 
by the TOTEM results at 7 TeV, Table V suggests that the case of fixed e = 0.080 presents better consistence 
with the experimental data than the case of fixed nig = 400 MeV. Since we shall not treat the simultaneous 
effects of both parameters, we may conclude that by fixing e at the above value we arrive at a best solution. 
Indeed this conclusion is corroborated by the following facts: 

• The bands for fixed nig lies, in general, inside those for fixed e and therefore the former may be seen 
as a particular case of the latter. 

• The fixed value e = 0.080 for the soft Pomeron intercept is consistent with all phenomenological analyses 
of the elastic hadron scattering, in particular with those treating constraint and extrema bounds [44, 45| . 

• The interval nig (MeV): [300, 600] is consistent with all phenomenological and theoretical results on 
the dynamical gluon mass. 

Therefore, we understand that our approach, as it has been presented here, has an optimum solution with 
only one fixed parameter, e = 0.080 and a relevant physical interval for the dynamical gluon mass scale, nig 
(MeV): [300, 600]. In this case, within the inferred error bands the predictions present reasonable consistency 
with the TOTEM results at 7 TeV, at least at the same level of some representative phenomenological models. 

4. Conclusions and Final remarks 

Further developments on a dynamical gluon mass approach to elastic and pp scattering have been 
presented and discussed. The main conceptual ingredient, introduced in 8, 9], concerns the dynamical gluon 
mass scale as a natural regulator for the IR divergences associated with the gluon-gluon cross section, leading 
to the identification of explicit nonperturbative contributions to high-energy elastic hadron scattering. In 
this context a consistent physical meaning is given to two ad hoc parameters, typical of mini-jet or QCD- 
inspired models: the infrared mass scale and the effective value of the running coupling constant. 

In this paper the main novel results concern the account of the energy dependencies in the dynamical 
gluon mass, Eq. (8), the detailed investigation on the influence in the predictions from relevant physical 
intervals for the dynamical gluon mass scale and the soft Pomeron intercept, as well as the introduction 
of a method to estimate the corresponding bands of uncertainty in all evaluated quantities. Fits to an 
improved data ensemble on pp and pp elastic scattering up to 1.8 TeV have led to quite good descriptions of 
all reduced data, including the high-energy differential cross section up to ^ 1.5 GeV^. With the inferred 
uncertainty bands the experimental data recently obtained by the TOTEM Collaboration are reasonably 
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well described, except the region beyond the dip position. Comparison with the predictions from some 
representative phenomenological models, presented and quoted in [1], shows that our results are at the same 
level of consistency with the TOTEM data, with the advantage in our case of explicit connections with 
nonperturbative QCD. As commented in Appendix A, our data reductions disfavor a soft solution for the 
dynamical gluon mass. 

The uncertainty in the evaluated quantities, associated with relevant intervals for the fundamental param- 
eters e and nig, seems to us a fundamental information for reliable phenomenological predictions. Although 
we did not address the question associated with both intervals simultaneously, we understand that our pro- 
posed estimation of the bounds gives more information on the reliability of the results than what could be 
obtained by fixing both parameters. Based on the results and discussion displayed in Subsection 3.3, we 
propose as our best present solution the case of fixed e = 0.080 and interval nig (MeV): [300, 600]. With 
that we have only one fixed parameter, with well founded physical justification. 

The critical point raised in respect to fixing parameters have been directed to a class of QCD inspired 
models. However the criticism certainly apply to any phenomenological model with fixed parameters whose 
numerical values do not have an explicit physical justification and its consequences in the evaluated quantities 
are not investigated or even discussed. 

In the phenomenological context, model developments are in general more important than any particular 
formulation and, in our case, at least two fundamental aspects demand further investigation. The first one 
concerns the phenomenological choice for the gluon distribution function, Eq. (18), which introduces the 
Pomeron intercept as a fundamental parameter. Despite its efficiency in the analysis we have developed, it 
may be important to study the effects of different parton density functions, along the lines that have been 
discussed by Achilli et al. [43]. In this class of QCD-based modelSjthe attenuation in the rise of the total cross 
sections comes from soft gluon emission from colliding partons [HI, @, 0| ■ These emissions affect the matter 
distribution and the resulting overlap function can be calculated from the soft gluon transverse-momentum 
resummed distribution. In this fc^ resummation scheme the overlap function in b-space depends upon the 
behaviour of the coupling as{kT) as fc^ — > 0, being the calculation only possible after the introduction of 
a phenomenological infrared modification for as{kT)- In the context of the DGM approach, it should be 
stressed that the QCD effective charge (7) is obtained from the QCD Lagrangian, i.e., it is derived from 
first principles, in a scenario where infrared effects are naturally taken into account. Hence this effective 
strong coupling could be naturally embodied in a soft fcy resummation scheme in order to avoid infrared 
divergences. A second aspect is related to the choice of form factors as simple dipole parametrizations. 
As commented before, the fact that our prediction for the exponent in the power law at large momentum 
transfer is near 10 (as is the case with the Aspen model), suggest that this result is directly related with the 
above choice for the form factors. Since deviations from the dipole parametrization have been indicated in 
model-independent analyses of pp elastic scattering (sil . Isil . [53| , to take into account these empirical results 
may give new insights in the predictions at large momentum transfer. We are presently investigating the 
above two lines. 

Appendix A 

With the DGM approach described in Sect. 2, in addition to the Cornwall's solution (8), we have also 
considered the possibility of a soft behavior for the dynamical gluon mass. A power-law running behavior 
for Mg{s) was first envisaged in [H^l and according to an OPE calculation the most probable asymptotic 
behavior of the running gluon mass is proportional to 1/s (HH ] . At the level of a non-linear Schwinger-Dyson 
equation this asymptotic behavior is given by 
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where 

4 6ci 

ci g [0.7, 1.3] and the parameters t and nig are constrained to he m the interval [1.0, 8.0] and [300, 600 
MeV], respectively 'H^. However, with the procedure described in Subsect. 3. A, this solution leads to much 
worse x^/DOF compared to the ones related to Eq. (8), which disfavor such solution. 
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